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Abstract 

CN 

A Fermi ball is a kind of non-topological soliton, which is thought to arise from the spontaneous 

CO 

breaking of an approximate Z2 symmetry and to contribute to cold dark matter. We consider a 



simple model in which fermion fields with multi-flavors are coupled to a scalar field through Yukawa 



Oh. coupling, and examine how the number of the fermion flavors affects the stability of the Fermi ball 



against the fragmentation. (l)We find that the Fermi ball is stable against the fragmentation in 
most cases even in the lowest order thin-wall approximation. (2)We then find that in the other 



> 
• 1— 1 

X 

specific cases, the stability is marginal in the lowest order thin-wall approximation, and the next- 
to-leading order correction determines the stable region of the coupling constants; We examine the 
simplest case where the total fermion number iV, and the Yukawa coupling constant Gi of each 
flavor i are common to the flavor, and find that the Fermi ball is stable in the limited region of the 
parameters and has the broader region for the larger number of the flavors. 

PACS numbers: 05.45.Yv, 95.35. +d 
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I. INTRODUCTION 

A Fermi ball 

1,2] 

, a kind of non-topological soliton [3(, is composed of three parts: a 
false vacuum domain, a domain wall enveloping the domain, and zero-mode fermions ^| 
confined in the domain wall. The Fermi ball is stabilized owing to the dynamical balance 
of the shrinking force due to the surface energy and the volume energy, and the expanding 
force due to the Fermi energy. The Fermi ball is thought to be a candidate for one kind of 
cold dark matter in the present universe j^|. 

Macpherson and Campbell pointed out that such stability holds good only for the spheri- 
cal shape of the Fermi ball [1]. They further showed that the Fermi ball is not stable against 
the deformation of the spherical shape, and thus flattens and fragments into tiny Fermi 
balls. The destabilization is caused by the volume energy of the Fermi ball. 

We, however, pointed out that the perturbative correction due to the domain wall curva- 
ture can stabilize the Fermi ball when the volume energy is small enough compared to the 
curvature effect 0, 3| . In case of a simple model with a single fermion flavor, we found that 
only in the quite narrow region of the parameters does the Fermi ball become stable. 

The purpose of the present paper is to examine how the fermion content of the model 
affects the stability of the Fermi ball. As an example, we consider an extended model in 
which fermions with multi-flavors are coupled to a scalar field through Yukawa coupling. 
Since the Pauli's exclusion principle does not apply to the different flavors of the fermions, 
the stable region of the parameters is expected to broaden. 



II. STABILITY OF FERMI BALL 



We consider the following Lagrangian density, 



1 n 

£ = ~(d,<P) 2 + E Wfy - GifiVi - U(<f>) , (1) 
A i=i 

where the scalar potential U(<f>) is given by 

W) = ^ 2 ~^) 2 + A(0) • (2) 

If the quantity |A(t>) — A(— v)\ is zero, the Lagrangian density is invariant under the Z2 
transformation, <-> —0. There is, however, a small but a finite quantity |A(u) — A(—v)\ ~ 
A <C Xv 4 , where the invariance is not a strict one. 



We consider a spherical Fermi ball with the radius R, and assume that the wave function 
and the boson <fi are static and that (f) depends only on the radial coordinate r. Let be 
the eigenfunction of the total angular momentum squared J , the z component J z and the 
parity P with the eigenvalues of J(J + 1), M and (— 1) J_< ^/ 2 (oj = ±1), respectively. Then, 
is written as 



' 1 g(r)y^(e^) 



(3) 



where yf§ and yjfj = {pxjr)yfj are the spherical spinors having the eigenvalues J and M, 
with J = I + u/2 = V — uj/2. Substituting Eq. (JHJ) into the Lagrangian L = J d 3 x C, we 
obtain 



dr 



where 



1 d 



K 



H f = °"i _ T~ + °"2 h v^Gi 

i dr r 



(4) 



(5) 



with = lo(J + |) and ipi(r) = (tyT)- Since the Fermi ball is a ground state with a fixed 
number of fermions, 

(6) 



N; 



d d x , 



we obtain the wave function ^ and the scalar field by extremizing 

L € [<f>, a] = m v<] + E e * IE /°° dr v-M - ^1 



(7) 



with the Lagrange multipliers The energy of the Fermi ball is expressed in terms of the 
fields as 



E 



dr 



Aitr^ 



2 V dr 



i KM 



(8) 



where is equal to the Fermi energy = dr ybjHfipi and is normalized as / °° dr ^1 V'i = 
1. In order to estimate the energy of the Fermi ball, we take the thin- wall approximation 
and obtain the correction due to the finite curvature radius R by the perturbation with 
respect to 1/R. We expand 0, ipi, and Hf in the power of 1/R, 



4> = 0o + 4>i H — 

^» = V'iO + H 

-H7 = i?Q + Hi + + ' ' • 



(9) 



where 



Ho = 0-1-3^- + &2^- + o- 3 Gi0o 
1 dr R 



H x 
H 2 



Kw 
Kw 2 



a 2 - 



R 3 



(10) 



with w = r — R. From 5L e /5<p = 8L e /8ipl = 0, we obtain the equations of motion, 



H i>io 



d 2 c 



dw 2 



+ E 



cj>=<t>0 



AnR 2 



E ^iO^iO 



(11) 
(12) 



K"Af 



and 



(H - e^ipa = -(H 1 - ea)ip. 



iQ 



' d 2 d 2 U 




2 d0 o 


dw 2 dip 2 


<P=4>o_ 


01= Rdw 



G 



E E Vta^i 



(13) 
(14) 



Neglecting A(0) in the scalar potential £7(0) for simplicity, we have analytic solutions for 

0o and Vio, 



4>q(w) = v tanh 
^io(w) 



1 1 



V^cosh^g X+ 



(15) 
(16) 

where 8^ = 2\~^v~ 1 is the thickness of the domain wall, 7$ = 2\~^Gi is the constant ^], 
Mi = dw cosh~ 27i u>/5& is the normalization factor, and x± 1S the eigenspinor of a 2 with 
the eigenvalue ±1. We note that the second term of the r.h.s. of Eq.(|12j) vanishes. The 
leading order of the eigenvalue is given by 

K 



ei0 ~R ' 

where we take K positive {u = +1). We have solutions for 0i and ipn, 

v w' „ h(w") 



(17) 



M w ) 



cosh 



—5 — / dw' cosh 4 — f dw'' 
h 2 ^ Jo 5 b Jo 



cosh 



2w h 



c l+ (w)x+ + Ci_{w)x- 



(18) 
(19) 
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where 



and 



[W) 6 b Rcosh 2 f b + 2nR^ ^ K L cosh 2 ^g 



cosh 7i f I FP Jo 5 b J w i cosh 2 ^ 







s b 

Substituting the solutions into Eq.(JSJ), we obtain the energy of the Fermi ball, 



C,_H = ^COsh--| dw'—^r 



E = Eq + 8E 



where E is the leading order contribution to the energy, 

SttAKAR 2 2£,A/. f 
-^0 = o r 



3 3R 

and (5-E 1 is the energy correction of the order of E Q x (5h/R) 2 , 



i 



T-N- 2 . f+°° w 2 



12R " J-oo cosh 4 f 

w' f w ' ^ „ h(w") 

b 



-2tt\iv 2 R / dw tt- / dw' cosh 4 — / dw'' 

J-oo COsh v- JO Oh Jo 



cosh 2 -^ 



3 

2 



2 „ iV/ /M-oo w 



3R 3 ^ j\fi J -oc cosh 27i f 



5 

,2 



5i? 5 V M ^~oo cosh 27l f io 6 b 
x / dw" 



x w" 



cosli 27 '^ 

+^E^cW dM) -T^ ( 



3i? 2Z r M 7-oo cosh 27! ^7o cosh 2 J 
o Ob Jo cosh ^7- 
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In the above equations, we use the relations, 



N i = E = E E=E (2^) = K max (K max + 1) , 

KM K=l M=-J X=l 



12 1 

^ t ^iO ~p) ^ t K K-max (K m ax 1 ) ( ^max q) 

KM ATM ^ 

2i\rl jV? 

(1) Stability in the leading order approximation 
Let us examine the stability of the Fermi ball within the leading order approximation in the 
£&/-R-expansion. From dE /dR = 0, we get the minimizing radius, 



Rmin l i j (^"0 

27T3 \6V 

and the energy at the radius, 

Eo = 2n*\vfeN^Yv . (28) 



we note d 2 E /dR 2 > at i? = -R mm - 

In order to examine the stability against the fragmentation, we compare two states; a 
state A in which a single Fermi ball has the fermion number Ni for z-th flavor, and a state B 
in which m Fermi balls have the fermion number each and conserve the total fermion 
number as J2™=i = Ni for each flavor. States A and B have the energy = E (Ni) 
and Eq = ^2 a Eo(N^ ), respectively. To compare the energy of the two states, we use 
Minkowski's inequality, 

(J2(NP + Nhif < (T(NP)if + fe(^ 2 ¥) S , (29) 

t i i 

where the equality is valid only for N- 2 ^ = cN^ (c > 0) with c being common for all i. 
Using the relation repeatedly, we have 

2 2 



EWJ <£(£(Af¥J , (so) 

where the r.h.s. is equal to the l.h.s. only for N^' = c^Ni (c^ > and Y^=i c ^ — !)• 
This leads us to the fact that except for the special case of = c^iVj, state A has lower 
energy than that of state B, and thus the Fermi ball is stable against the fragmentation in 
the leading order approximation. This situation that the Fermi ball is stable in most cases 
is characteristic of the case with multi-flavor of fermions, and qualitatively different from 
the case of a single flavor []]]. In case of = c^A^, the two states have the same energy 



in the leading order approximation, and the correction term 5E determines the stability of 
the Fermi ball against the fragmentation. 

(2) Stability in the next-to-leading order approximation in the special case 

N M = c (a) N . 

We examine the stability of the Fermi ball against the fragmentation in the case of iV"/ = 
c^Ni. Substituting R = R min into Eq.flMJ) yields 

5E = C(X,G i ,N i )v , (31) 

where 

L{A,Lri,l\i) — g— j 1 —j 1 5 

6(E,^)3 3\*(EiN?) 
2048vr§ &iN?WMi)) + 128tt(E, GjNfWMi)) 
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r, 



5A6(£^) 3 3A(E,iV/) 

Here, I\ to I5 are given by 



dx 



00 
-i- 00 



cosh 4 x 



h = I dx f dx' coshV f ' dx" 

i-00 cosh x J jo cosh x" 



00 J 



/ 3 (Y) = / dx 



00 cosh 27l x 

dx — —~- — / dx' cosh 27i x' / dx" — -5- — 
-00 cosh 7l x Jo Jx 1 cosh 7 *x" 

/ + OO rp fX i fX 1 

dx — \ dx ' — -2- / dx " coshV 
-00 cosh 7l x Jo cosh x' Jo 

x f X " dx"' h (?P , (33) 
Jo cosh x 

with rescaled as h(x) = ^ k h(5bx) and Mi as JVi = ^-M- AVe compare state ^4 of 

the single Fermi ball and state B of m Fermi balls with the total fermion number to be 
conserved for each flavor. States A and B have the energy Ej± = Eq(N{) + C(X,Gi, N^v 
and E B = £ a E (N^ a) ) + Ea C(A, G 4 , N^ a) )v, respectively. In case of N$ a) = c^i^, we derive 
E a E (N^ a) ) = E (Ni) from Eq.(j2HD and C(A, G i; ivf >) = C(X, G u N { ) from Eq.flH2J), and thus 
find that state i3 has the energy Eg = E (Ni) + mC(X,Gi, Ni)v. Therefore, if C(X,Gi, Ni) 
is positive, state A has lower energy than that of the state B by the magnitude of the 
correction term SE, and the Fermi ball is stable against fragmentation even in the special 
case of iV/ a) = c^iV,-. 
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Let us consider the simplified model to examine how the number of the fermion flavors 
n affects the stability of the Fermi ball in case of iVf } = c^Ni. We assume that belongs 
to a multiplet of the internal symmetry with a common Yukawa coupling constant G and 
also assume that the fermion number is common to the flavor, i.e., iVj = N. Under these 
assumptions, the coefficient C is independent of N and dependent on A, G and n from 
Eq.(jnU). We evaluate Eq. (JH2j) using a numerical integration, and obtain the stable region of 
the parameters where C is positive (see Figures^and|21). These figures show that the allowed 




FIG. 1: The allowed regions (shadowed) of the scalar self-coupling constant A and the Yukawa 
coupling constant G for the Fermi ball to be stable against the fragmentation. We assume that 
the fermion (1 < i < n) belongs to a multiplet and the boson eft to a singlet of the internal 
symmetry, and that the fermion number iVj is common to the flavor as N{ = N. The figure shows 
that the allowed region broadens as n increases. 

regions of the parameters exist for the Fermi ball to be stable against the fragmentation (the 
shadowed regions in the figures). We see in the figures that the allowed region broadens as 
the number of the flavors n increases. 
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FIG. 2: The allowed region (shadowed) of the scalar self-coupling constant A (left) and the Yukawa 
coupling constant G (right) for the Fermi ball to be stable. The assumptions are the same as those 
in Figure ^ We see that the allowed regions broaden as n increases. 

III. CONCLUSION 

We have considered a model for the Fermi ball in which the fermions with multi-flavors 
(1 < i < n) are coupled to the scalar field <fi and the total fermion number of each i-th 
flavor is fixed as iVj. We have examined the region of the parameters for the Fermi ball to 
be stable against fragmentation, and how the number of the fermion flavors n affects the 
stability. 

We have considered the thin- wall Fermi ball, i.e., the radius R is much larger than the 
wall thickness 5b- We have taken into account the effect due to the finite wall thickness by the 
perturbation expansion with respect to 5b/ R- In the leading order thin- wall approximation, 
we have compared the energy of the initial state of a single Fermi ball and that of the 
final state of fragmented m Fermi balls, with the total fermion number Ni of each flavor 
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i being conserved, J2™=i N- a * = N{. We have found that the former is smaller than the 
latter and thus the Fermi ball is stable against fragmentation, except for the special case 
of = c^Ni with Y^=\ c ^ — 1- This situation that the Fermi ball is stable in most 
cases is characteristic of the case with multi-flavor of fermions, and qualitatively different 
from the case of a single flavor. In the special case of = c^iVj, the two states have the 
same energy in the leading order approximation and the next-to-leading order correction 
term SE determines the stability. There we have found that the energy of the initial state 
is Eq + Cv and that of the fragmented states is E + mCv, where v is a symmetry breaking 
scale and C is a coefficient dependent on the scalar self-coupling constant A, the Yukawa 
coupling constant Gi and the fermion number N%. This tells us that even in that case the 
Fermi ball is stable when C takes a positive value in the parameter region of A, Gi and iVj. 

We have considered the simplified model in which a multiplet of fermions has a common 
Gi and a common iVj for each flavor i. We have found that the allowed region of the 
parameters for the Fermi ball to be stable exists and broadens as the multiplet dimension n 
increases. 
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